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Separated before birth: pulsars B2020+28 and B2021+51 as the 

remnants of runaway stars 
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Abstract. Astrometric data on the pulsars B2020+28 and B202 1+51 suggest that they originated within several parsecs of each 
other in the direction of the Cyg OB2 association. It was proposed that the pulsars share their origin in a common massive binary 
and were separated at the birth of the second pulsar following the asymmetric supernova explosion. We consider a different 
scenario for the origin of the pulsar pair based on a possibility that the pulsars were separated before their birth and that they are 
the remnants of runaway stars ejected (with velocities similar to those of the pulsars) from the core of Cyg 0B2 due to strong 
three- or four-body dynamical encounters. Our scenario does not require any asymmetry in supernova explosions. 
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1. Introduction 



Most stars form in dense embedded clusters and reside in bi- 
nary systems (either primordial or tidal). If both binary compo- 
nents are massive enough, they end their lives as core-collapsed 
supernovae (SNe). Stellar remnants of SN explosions, usually 
the neutron stars (NSs), have peculiar velocities at least an or- 
der of magnitude higher than the typical velocities of their pro- 
genitors, the OB stars (e.g. Gunn & Ostriker [1 9701 1 . It is be- 
lieved that the high velocities of NSs are due either to the asym- 
metry of SN explosions (e.g. Dewey & Cordes 198 7 ) or to the 
disruption of tight massive binaries following the second (sym- 
metric) SN explosion (e.g. Iben & Tutukov|T996|l. The progress 
in measuring the proper motions and parallaxes of NSs (pul- 
sars) allows their pecuUar (transverse) velocities to be deter- 
mined with high precision and makes it possible to trace their 
trajectories back to the parent star clusters (e.g. Hoogerwerf et 
al. 120011 Chatterjee et al. I2005I I. Recently VIemmings et al. 
( 12004 1 hereafter VCC04) have used the high-precision astro- 
metric data (proper motions and parallaxes) for two dozen pul- 
sars to determine their trajectories in the Galactic potential and 
to search for pairs with a common origin. They discovered that 

o 

two pulsars from their sample (presently separated by ~ 23 ) 
originated within several parsecs of each other in the direction 
of the Cyg OB2 association. VCC04 interpret their discovery 
as an indication that the progenitors of both pulsars, B2020+28 
and B2021+51, were the members of a common massive bi- 
nary and suggest that the pulsars were separated at the birth of 
the second one following the asymmetric SN explosion. 



In this Letter we explore a different scenario for the ori- 
gin of B2020-H28 and B202 1 -1-5 1 . We suggest that these pulsars 
were separated before their birth and that they are the remnants 
of runaway stars ejected (with velocities similar to those of the 
pulsars) from the parent star cluster due to the strong three- or 
four-body dynamical encounters. Our scenario does not require 
any asymmetry in SN explosions. 

2. Pulsars B2020+28 and B2021+51 : origin in a 
common binary 

The main result presented in VCC04 is that B2020+28 and 
B2021+51 originated within several parsecs of each other. 
VCC04 derived the most likely three-dimensional peculiar ve- 
locities of the pulsars at birth, ^150 and ^ 500 km s"' (respec- 
tively, for B2021+51 and B2020-H28), and the angle between 

o 

the velocity vectors i/r ^ 160 . These velocities can, in prin- 
ciple, be produced via disintegration of a tight (semi-detached) 
massive binary after the second (symmetric) SN explosion (e.g. 
Iben & Tutukov 1996 ); in this case, however, if/ is always < 90 . 
Moreover, the spin characteristics of B2020-I-28 and B202 1 H-5 1 
(typical of non-recycled pulsars) argue against the origin of 
these pulsars in a common tight binary (cf. VCC04). 

One possible way to reconcile the kinematic data with the 
common binary scenario is to assume that the binary was dis- 
rupted either after the first or the second asymmetric SN ex- 
plosion (VCC04). Note that the similarity between the pulsar's 
characteristic ages (^ 2.88 and ^ 2.75 Myr) impUes that the 
mass ratio of the binary components was ~ 1. Therefore, de- 
pending on the initial parameters (binary separation, etc), the 
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Fig. 1. The dependence of the velocities of the (first-born) pul- 
sar and its former companion star (now the runaway progenitor 
of the second pulsar) on the angle between the kick vector and 
the direction of motion of the exploding star (shown, respec- 
tively, by the solid and the long-dashed lines). The horizontal 
short-dashed lines indicate the pulsar velocities suggested by 
VCC04. See text for details. 

binary system at the moment of the first SN explosion consists 
of two red supergiant or Wolf-Rayet stars or of two carbon- 
oxygen (CO) cores. 

The latter situation can be realised if the massive binary 
evolves through two common-envelope phases (see Belczyriski 
& Kalogera 2001 ). A natural outcome of the evolution of this 
type of binaries, provided that the SN explosions were of zero 
or moderate asymmetry, is the origin of a binary non-recycled 
pulsar (Belczynski & Kalogera '2001V The CO binary, how- 
ever, could be disrupted after the first (or the second) asym- 
metric SN explosion if the kick received by the stellar remnant 
was of proper magnitude and orientation (see Tauris & Takens 
fT998l l. 

For illustrative purposes, we consider the disruption of a 
CO binary following the first asymmetric SN explosion. For 
parameters of the CO binary given in Belczyriski & Kalogera 
( l200TT l and using Eqs. (44)-(47) and (51)-(56) given in Tauris 
& Takens (1998), one can show that the pulsar velocities and i/^ 
could be explained if the kick imparted to the first-born pulsar 
(B2020-I-28) was ~ 500 km s"' (for the sake of simplicity we 
assume that the second SN explosion was symmetric), while 
the angle between the kick vector and the direction of mo- 
tion of the exploding star, 6, was ^ 4013 (see Figs. 1 and 2 
and Gvaramadze 12006 1. It is obvious that the kick should be 
stronger if, at the moment of the first SN explosion, the binary 
consists of red supergiant or Wolf-Rayet stars (cf. VCC04). 

Another possibility is that the pulsars attained their veloci- 
ties in the course of disintegration of the binary after the second 
asymmetric SN explosion. Since both pulsars are not recycled, 
one should assume either that the binary separation was suf- 
ficiently large (so that the wind of the secondary star did not 
affect the evolution of the first-born pulsar) or that the binary 
evolved through a double common-envelope phase (see above). 
VCC04 suggest that the pulsars were born in a wide binary, but 
in their analysis they draw an erroneous conclusion that the pul- 

' Note that for 64 i 6 i 290 , the binary system remains bound. 




Fig. 2. The angle between the velocity vectors of the first- and 
second-born pulsars as a function of the angle between the kick 
vector and the direction of motion of the exploding star. The 
horizontal dashed line indicates the angle between the pulsar 
velocity vectors suggested by VCC04. 

sar velocities can be explained by a kick of only ^ 200 km s ' 
(see Gvaramadze I2006l l. One can show, however, that in both 
the above-mentioned cases the kick imparted by the second SN 
explosion should be > 500kms 

Thus, the origin of the pulsars in a common binary implies 
that at least one of the SN explosions was asymmetric enough 
to produce a kick of > 500 km s ' . If, however, SNe can indeed 
impart high velocities to NSs, then it is not necessary to assume 
that the pulsars originated in the same binary, but instead one 
can suggest that they were created by two separate SN explo- 
sions occurred in the same parent star cluster within a few 10^ 
yr. Our scenario for the origin of the pulsar pair has something 
in common with the latter possibility, but we do not require any 
asymmetry in the SN explosions. 

3. Pulsars B2020+28 and B2021+51 : dynamical 
ejection from the young massive star cluster 

The recent discovery of the so-called hypervelocity stars 
(Brown et al. 2005) and hyperfast pulsars (Chatterjee et al. 
.2005) . the ordinary stars and pulsars moving with extremely 
high (~ 1 000 km s"') peculiar velocities, suggests a possibility 
that the hypervelocity stars could be the progenitors of hyper- 
fast NSs, provided that they are massive enough (Gvaramadze 
et al. I2007I I. A strong argument in support of this possibility 
comes from the fact that the mass of one of the hypervelocity 
stars, HE 0437-5439, is ^ SMq (Edelmannet al. .2005j so that, 
in principle, it can end its evolution as a hyperfast NS ! The high 
velocities (~ 200 - 400 km s"') inferred for some early B-type 
stars at high galactic latitudes (Ramspeck et al.|2001) also sup- 
port the possibility that high-velocity pulsars could originate 
from high- velocity runaway stars. 

Gvaramadze et al. (2007) suggest that the origin of hy- 
pervelocity stars could be connected not only with scatter- 
ing processes involving the supermassive black hole (BH) in 
the Galactic centre (the common wisdom; originally suggested 
by Hills [19881 see also Yu & Tremaine l2003l Gualandris et 
al. '2005), but also with strong three- or four-body dynami- 
cal encounters in the dense cores of young massive star clus- 
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ters (YMSCs), located either in the Galactic disk or near the 
Galactic centre. The discovery of a halo population of early B 
stars, whose lifetimes are shorter than the times-of-flight from 
the Galactic centre (Brown et al. ,2007. see also Ramspeck et 
al. 1200 IK supports this suggestion. We believe, therefore, that 
the pulsars B2020+28 and B2021+51 could be the remnants of 
high-velocity runaway stars ejected from the same YMSC. The 
kinematic and characteristic ages of the pulsars (respectively, 
~ 2 and ~ 3 Myr; VCC04) imply that by the moment of ejec- 
tion the progenitor stars have already become Wolf-Rayet stars 
[the short-lived (< 1 Myr) helium (He) cores of massive stars; 
cf . Gvaramadze et al. 120071 . 

Of the two mechanisms that could be responsible for 
the origin of the high-velocity progenitors of B2020H-28 and 
B202 1-1-51, the first reUes on close dynamical encounters be- 
tween hard (Heggie 1 19751 1 massive binary stars in the dense 
core of a YMSC. The peculiar velocities of runaway stars pro- 
duced in this process are similar to the orbital velocities of the 
binary components (e.g. Leonard & Duncan ri990l l. but occa- 
sionally they could be much higher. Scattering experiments by 
Leonard ( I1991I I showed that the maximum velocity attained by 
the lightest member of the binaries involved in the interaction 
(e.g. the He core of a massive star or an early B-type star) can 
be as high as the escape velocity, Vesc, from the surface of the 
most massive star in the binaries. 

For the main-sequence stars with the mass-radius relation- 
ship (Habets & Heintze [MB, rus = 0.8(mMs/Mo)°-^ /?o, 
where tms and niMS are the stellar radius and the mass, the 
maximum possible velocity of ejected stars is a weak function 
of otms, V™" ^ Vesc - 7OOkms-'(mMs/M0)" '5 and could be 
as high as ~ 1 400 kms"' (cf. Leonard 1991). Numerical sim- 
ulations by Leonard (.1991^ showed that about 4% of binary- 
binary encounters result in the origin of runaway stars with 
Vej - 0.5 Vesc, which is enough to explain the velocity of 
~ 500 km s"' suggested by VCC04 for one of the pulsars. Note 
that the results of Leonard ( |1991| l were used by Tenjes et al. 
( 1200 11 1 to explain the origin of the high-velocity (~ 400 kms') 
runaway star HIP 60350. 

Another possible mechanism for producing high-velocity 
stars is based on exchange encounters between tight binary 
stars and a compact massive object, either a very massive 
star (VMS), formed through the runaway stellar collisions and 
mergers in the mass segregated core of a YMSC (e.g. Portegies 
Zwart et al. '1999), or its descendant, an intermediate-mass BH 
(e.g. Portegies Zwart & McMillan I2002I I. After the close en- 
counter and tidal breakup of the binary, one of the binary com- 
ponents (usually the more massive one) becomes bound to the 
compact object, while the second one recoils with a high veloc- 
ity given by Vej ~ [M/(mi -H OT2) ]'^''(2GOTi/a)'^^ (Hills 119881 
see also Gvaramadze et al. 120071 ). where M is the mass of the 
compact object, nii and ma are the masses of the binary com- 
ponents (mi > ma), and a the binary separation. 

In YMSCs of mass ~ IO'^Mq, the mass of the VMS 
does not exceed several 100 M©, while the thermal (Kelvin- 
Helmholtz) time scale of the VMS is shorter than the mean 
time between collisions (see Portegies Zwart et al. I1999I I. In 
this case, the growing VMS rapidly evolves to the thermal equi- 



librium (e.g. Suzuki et al. '20071), so that one can adopt the fol- 
lowing mass-radius relationship, R ^ 1.6 (M/Mq)*' ''^TJq, where 
R is the radius of the VMS (see Freitag et al. 120061 and refer- 
ences therein). In the process of an exchange encounter with 
a binary, the VMS could be considered as a point mass if the 
binary tidal radius, ~ [M/(mi + m2)]'^"'fl, is at least several 
times larger than R. For M = 200 - 300Mo,mi = 30 Mq (a 
main-sequence star), m2 = 8 Mq (a He core), and a = 50 Rq, 
one has rt 90 - 100 (i e. much larger than ^ 19-23/?0) 
and Vej ^ 630 -670 kms"', that is enough to explain the pulsar 
velocities. 

In more massive (> 10'' Mq) YMSCs, the VMS can acquire 
a mass of several 1 000 Mq (Portegies Zwart et al. 120041 ). But 
in this case, the thermal time scale is comparable to the colli- 
sion time and the VMS remains inflated untill collapsing into 
an intermediate-mass BH (e.g. Portegies Zwart et al. 120061 ). 
Exchange encounters with this VMS would not produce high 
ejection velocities. The star ejection from the YMSC, however, 
would be very effective if the VMS leave behind a BH of mass 
~ 1 000 Mq (e.g. Gualandris & Portegies Zwart 12007] ). 

4. Cyg 0B2 

The astrometric data on B2020+28 and B202 1-1-51 suggest that 
these pulsars (or their progenitors; our preference) were ejected 
~ 1.9 Myr ago from the same origin at a distance of ~ 1.9 
kpc in the direction of the Cyg OB2 association (VCC04). 
The parent YMSC (or its descendant) should still be located 
at about the same distance since its possible peculiar velocity 
of < 30 kms ' (typical of the OB associations near the Sun; de 
Zeeuw et al. 1999) would result only in a slight offset of < 60 
pc (cf. VCC04). To constrain the current age of the parent clus- 
ter, we assume that the initial mass of the progenitor stars of 
B2020H-28 and B2021+51 could be as high as > 50 Mq. (It is 
believed that stars of this mass can lose most of their mass via 
stellar wind or mass transfer on a binary companion and leave 
behind NSs; e.g. Vanbeveren et al. ,1998; Wellstein & Langer 
119991 cf. Muno et al. 120061 ) From this it follows that die mini- 
mum age of the parent YMSC should be ~ 5 Myr, that is, ~ 2 
Myr (the pulsar kinematic age) plus ~ 3 Myr (the lifetime of a 
^ 50 Mq star). Assuming that the YMSC initially contained at 
least 10 stars of mass > 5OM0, one has the (initial) mass of the 
cluster of > 10"* Mq (for a 0.2 - 120 Mq Salpeter initial mass 
function). 

The only likely candidate for the birth cluster of B2020-I-28 
and B202 1-1-51 in the region suggested by VCC04 is the Cyg 
OB2 association. Numerous star clusters in its neighbourhood 
(see, e.g., Le Duigou & Knodlseder '2002) cannot pretend to 
play this role either due to their youth or low masses. 

Cyg OB2 is one of the most massive and compact OB as- 
sociations in our Galaxy (Knodlseder 2000). The large number 
(~ 100) of O stars identified in Cyg OB2 (Knodlseder 2000, see 
also Comeron et al. 120021 ) implies that its mass could be as high 
as ~ 10 Mq. The angular radius of Cyg OB2 is ~ 1 , while the 
half light radius is ~ 13 (Knodlseder .2000) , that at the dis- 
tance of Cyg OB2 of ~ 1.5-1.7 kpc (Hanson |2003l Massey 
& Thompson 1991) coiTesponds, respectively, to ~ 25 - 30 pc 
and ~ 5 - 6 pc. Note that the centre of Cyg OB2 lies within the 
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2cr likelihood contour of the pulsar birth location and, at the 2cr 
level, the distances to the Cyg 0B2 and the birth location are 
consistent with each other Age estimates for Cyg OB2 range 
from ~ 1 to 5 Myr (e.g. Bochkarev & Sitnik |1985l Knodlseder 
et al. |2002| l. The wide age spread suggests that the star forma- 
tion in the Cyg OB2 was non-coeval. The non-coevality could 
be understood if the star formation in the association started 
initially in the dense core of the parent molecular cloud and 
then propagated to its lower density periphery. It is believed 
(e.g. Elmegreen l2000l ) that the star formation occurs on one or 
two dynamical time scales, fdyn ~ (Gp)"'^^, where p is the gas 
density in the cloud, so that in a density-stratified cloud of mass 
of ~ 10 times higher than the stellar mass of Cyg OB2 and the 
size similar to that of the association, the age spread could be 
comparable with fdyn ~ 5 Myr. 

If the progenitor stars of B2020-I-28 and B202 1-1-51 were 
ejected from Cyg OB2, then we suppose that, ~ 2 Myr ago 
(or ~ 3 Myr after the formation of the first massive stars in 
the centre of the association), the core of the association was 
much more compact and denser. Assuming that the associa- 
tion expands with a velocity equal to its velocity dispersion 
(^ 2.4kms Kiminki et al. I2006I I. one finds that the stars 
located within the current half light radius were originally con- 
centrated in a region of a radius of < 1 pc. It is likely that the 
two star clusters projected close to each other (~ 6 ) near the 
centre of Cyg OB2 (Bica et al. 120031) are the remainders of 
this dense core. We suggest that the mass of the core was much 
higher (> lO'* Mq) than the current mass of the clusters (several 
1 OOOMq; Bica et al. 120031) and that it was significantly reduced 
during the last 2 Myr due to the overall expansion of the associ- 
ation and star ejections following close dynamical encounters 
[recent N-body simulations by Pflamm-Altenburg & Kroupa 
(1200 6') showed that dynamical processes in the Trapezium clus- 
ter could be responsible for the loss of at least 75% of its initial 
content of OB stars]. Thus we believe that ~ 2 Myr ago the 
conditions in the core of Cyg OB2 were favourable to the dy- 
namical processes discussed in Sect. 3. 
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